function function_f=Phi_fun(variables_v,parameters)
BETTA=parameters(1);
B=parameters(2);
CHI=parameters(3);
CHI0=parameters(4);
THETA=parameters(5);
DELTA=parameters(6);
ALFA=parameters(7);
PHI=parameters(8);
PHIzero=parameters(9);
ZI=parameters(10);
KAPAw=parameters(11);
ETA=parameters(12);
RHOR=parameters(13);
PHIpai=parameters(14);
PHIpai_1=parameters(15);
PHIy=parameters(16);
PHIy_1=parameters(17);
PHIc=parameters(18);
PHIc_1=parameters(19);
PHIl=parameters(20);
NU=parameters(21);
U0=parameters(22);
U0d=parameters(23);
RHOz=parameters(24);
RHOd=parameters(25);
RHOn=parameters(26);
RHOp=parameters(27);
RHOa=parameters(28);
Kss=parameters(29);
OUTPUTss=parameters(30);
PAIss=parameters(31);
MUZss=parameters(32);
Css=parameters(33);
AA=parameters(34);
lss=parameters(35);
Rss=parameters(36);
Wss=parameters(37);
c_ba1=variables_v(1,:);
muz_cu=variables_v(2,:);
d_cu=variables_v(3,:);
n_cu=variables_v(4,:);
paistar_cu=variables_v(5,:);
a_cu=variables_v(6,:);
sigma_perturbation=variables_v(7,:);
function_f=[];
